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NOMENCLATURE 


ACA 

Attitude Controller Assembly 

ACB (X, Y, Z) 

Accelerometer bias (channels X, Y, Z) 

ADSRA (X, Y, Z) 

Gyro drift due to acceleration along the 
sp!n reference axis (Channels X, Y, Z) 

ADIA (X, Y, Z) 

Gyro drift du® to acceleration along the 
input axis (Channels X, Y, Z) 

AGS 

Abort Guidance System 

AOT 

Alignment Optical Telescope 

APS 

Ascent Propulsion System 

ASA 

Abort Sensor Assembly 

BDA 

Bermuda (tracking station) 

CDU 

Coupling Data Unit 

COAS 

Crew Optical Alignment Sight 

CM 

Command Module 

CMC 

Conmand Module Computer 

DAP 

Digital Auto Pilot 

DEPA 

Data Entry and Display Assembly 

DPS 

Descent Propulsion System 

DSKY 

Display and Keyboard 

EPC 

Earth Prelaunch Calibration 

FDAI 

Flight Director Attitude Indicator 

GDA 

Glmbal Drive Actuator 

GDS 

Goldstone (tracking station) 

GET 

Ground Elapsed Time (Range Time) 

G&N 

Guidance and Navigation 

GSOP 

Guidance System Operational Plan 

HOPE 

Houston Operations Predictor/Estimator 

IFC 

Inflight Calibration 

IMU 

Inertial Measurement Unit 

TRIG 

Inertial Rate Integrating Gyro 

JPl 

Jet Propulsion Laboratory 
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LGC 

LM Guidance Computer 

LM 

Lunar Module 

LOS 

Line-of-sight 

LOX 

Liquid Oxygen 

MCC 

Midcourse Correction 

MERU 

Mi lli -Earth Rotational Units 

MIC 

Minimum Impulse Control (mode) 

MSC 

Manned Spacecraft Center 

NASA 

National Aeronautics and Space Administration 

NAT 

NASA Apollo Trajectory 

Omega P' error 

Rate error about P axis 

Omega U' error 

Rate error about U‘ axis 

Omega V' error 

Rate error about V' axis 

P error 

Yaw axis error 

PGNCS 

Primary Guidance, Navigation and Control 
System 

PIC 

Pre-Installation Calibration 

PIPA 

Pulsed, Integrating Pendulous Accelerometer 

PPM 

Parts per Million 

PTC 

Passive Thermal Control 

RCS 

Reaction Control System 

REFSMMAT 

Reference to Stable Member Matrix 

RHC 

Rotational Hand Controller 

RSS 

Root of the Sum of the Squares 

RTCC 

Real Time Control Center 

S/C 

Spacecraft 

SFE ( X , Y, 2) 

Scale Factor (Channels X, Y, Z) 

SM 

Service Module 

SODB 

System Operational Data Book 

SPS 

Service Propulsion System 

SIVBIU 

Saturn IVB Instrumentation Unit 
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Nomenclature (Continued) 


TEI 
THC 
TLI 
TTCA 
U error 
U' error 
V error 
V' error 
VG 

VO (X, Y, Z) 


Trans-Earth Insertion 

Translational Hand Controller 

Trans-Lunar Injection 

Thrust and Translation Controller 

Computed Errors 

Computed Errors 

Computed Errors 

Computed Errors 

Velocity Gained 

Velocity Offset (X, Y, Z) 

Measured gravity vector in IMU coor- 
dinates (x) 

Measured gravity vector in IMU coor- 
dinates (y) 

Measured gravity vector in IMU coor- 
dinates (z) 

Micro-gravities 


P, U, V axis DAP control axis oriented relative to 

LM body axes as shown below: 
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1.0 INTRODUCTION 


This report presents the conclusions of the analyses of the inflight 
performance of the Apollo 13 mission Guidance, Navigation and Control 
equipment onboard the CSM-109 and IM-7 spacecrafts. This analysis is 
supplanent 1 to the Apollo 13 Mission Report (Reference 1) . 
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2.0 SUMMARY 


CM IMU 

Analysis of the ascent and TLI bum errors indicated an X-gyro con- 
stant drift which was outside of the la specification (2 meru). The 
derived values were 3.14 meru for ascent and 2.35 meru for TLI. There 
is evidence that this apparent instrument error actually resulted from 
an inappropriate gyro drift compensation load. PIPA bias values were 

reasonably stable during the time the IMU was turned on, but the Z PIPA 

2 

bias value shifted approximately 1.64 cm/sec across the long power down 
period and required a CMC compensation update prior to entry. 

After the SM LOX tank incident occurred, the IMU power supply under- 
went degradation. Concurrently, the PIPA's registered a 1c level 
acceleration which was first interpreted as venting. Howeve. , doppler 
radar failed to confirm venting of the necessary magnitude. It is now 
believed that the degradation in power supply voltage induced a corres- 
ponding transient in PIPA performance. 

CM OPTICAL NAVIGATION 

Processing of translunar P23 (star-horizon) data indicates that the 
actual horizon altitude was about 18 kilometers. However, the apparent 
altitude (due to small sighting errors) was about 10 kilometers. 

LM IMU 

IMU performance was good. PIPA bias values were quite stable about 
their prelaunch values. No direct measurement of gyro drift was obtained, 
but there is good evidence that total (RSS) misalignment of the IMU at 
the time of the TEI bum-after 20.5 hours of gyro drift - was of the 
order of 0.5 to 0.6 degrees. The lo value from LM IMU drift alone 
(ignoring Initialization errors) is 1.1 degrees, indicating excellent 
gyro performance. 



LM DIGITAL AUTOPILOT 


During the mission, the LM DAP was called upon to control both LM/CSM 
and LM/CM spacecraft configurations in auto maneuver and attitude hold modes. 
Performance was satisfactory in all cases although in some instances 
violations of deadband limits did occur. 

The roll GDA was observed to drive approximately -1.3 degrees from 
its initial position at the start of the TEI burn, stimulating conjecture 
that the engine gimbal trim function might have been abnormal. Detailed 
analysis indicated that this position change was necessary to relieve com- 
pliance and correct for initial mistrim. Performance is now believed to 
have been nominal. 

Difficulty was encountered when attempting to maneuver the space- 
craft into the PTC attitu^® following the TEI burn. Downlink data were 
not available for the maneuver, so that investigation was necessarily 
restricted. One explanation was found by a theoretical examination of cross 
couplings resulting from various control modes. It was determined that 
the use of rotational commands would have resulted in significant cross 
coupling (due to jet impingement forces) and drastic alterations in the 
intended commands. Although purely translational commands were planned for 
this maneuver, it is hypothesized that some rotational commands were inter- 
mingled and were the source of the problem. 

Another possible explanation is that the difficulty arose from the 
necessity of determining hand controller conmands purely by interpreting 
gimbal angle displays on the DSKY. This was necessary because the FDAI was 
powered down. 

LM ABORT GUIDANCE SYSTEM 

The AGS was used for spacecraft control In two burns, MCC 3 and 
MCC 4. In both cases, performance was satisfactory. 
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AGS gyro and accelerometer errors were estimated from free flight 
data and from AGS/PGNCS velocity comparisons obtained during the TEI 
burn. Instrument static errors showed excellent stability. The TEI 
burn afforded the only opportunity for observing acceleration sensitive 
errors. Instrument dynamic performance during that burn was within the 
2 o limits determined from the ASA023 error model. This is excellent per- 
formance, particularly in view of the fact that the ASA023 dropped 

o 

23 degrees F below its specified minimum temperature (60 F) during the 

24 hour period in which it was shut down prior to the TEI burn. 

LM OPTICAL ALIGNMENT CHECKS 

Prior to the TEI burn the LGC was used to aim the AOT line-of-sight 

in the LGC calculated sun direction. This was done to assess IMU align- 
ment errors and resulted in an estimated misalignment of 0.5 degrees. 

By its nature, the check could not resolve errors about the line-of-sight. 
However, the line-of-sight direction was such that misalignments about 
it had negligible impact on significant state vector errors induced by 
the IMU during the TEI burn. As conducted, the check determined only 
the magnitude of the misalignment in one plane - that plane in which 
misalignments were of greatest consequence. The direction of these 
errors was not determined, and therefore the appropriate corrective IMU 
torques could not be calculated. The IMU misalignment generated a burn 
error of approximately 5.3 feet/second in the TEI burn. Unfortunately, much 
of this error was in the direction most critical to entry flight path angle; 
the resultant (potential) error in the post-TEI trajectory entry flight 
path angle was in excess of -4 degrees. The allowable error in flight path 
angle is 0.5 degrees. Two midcourse corrections were subsequently required 
to correct this error. 
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An alignment of the LM IMU was performed using the sun and moon 
as optical targets in preparation for the MCC 4 burn. A star angle dif- 
ference of -1.1? degrees was computed by the LGC and displayed on the 
DSKY, indicating a very large astronaut sighting error. It was sub- 
sequently determined that the actual sighting error was only about .08 
degree. The remainder of the apparent sighting error arose from in- 
accurate LGC ephemeris data for the moon and sun, and from the LGC 
software. 
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3.0 CM SYSTEMS 


3.1 CM INERTIAL MEASUREMENT UNIT 


3.1.1 Velocity Comparisons During Ascent and TLI 


The Apollo 13 CM IMU performance analysis was based on comparisons 
of Apollo (denoted G&N) and Saturn (denoted S-IVB) measured velocities. 
Analysis centered around the "sensed" velocities - those resulting 
from integration of that portion of the vehicles acceleration which can 
be sensed by the accelerometers. Sensed velocities exclude the influence 
of the gravitational field. In addition to the sensed velocity com- 
parison, a cross check of ascent phase results was obtained by analyzing 
"total" velocity differences. Total velocity is the actual vehicle 
velocity in inertial space. It is obtained by integrating both sensed 
acceleration (due to engine thrusting) and gravitational acceleration. 

The cross check was performed because telemetry data dropouts caused 
minor discrepancies in the GN sensed velocity estimates. Sensed and 
total velocity differences are presented below for the end of the ascent 
phase (t = 752 seconds, CMC clock time) and for TLI cutoff (t = 9704.48) 

Type of aX aY aZ 

Time Comparison Comparison (ft/sec) (ft/sec) (ft/sec) 


752.0 sec 
752.0 sec 
9704.48 sec 


Sensed 

G&N - S-IVBIU 

-6.75 

87.04 

- 3.82 

Total 

G&N - S-IVBIU 

-5.41 

73.51 

1.02 

Sensed 

G&N - SIVBIU 

-20.50 

- 2.50 

2.03 


IU represents the edited Saturn telemetry data. 
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The ascent and TLI uncompensated sensed velocity differences (G&N - 
S-IVBIU) appear in Figures 3-1 through 3-6. Those IMU errors solved for 
in this analysis are defined in Table 3.1. IMU error sets derived to 
fit the ascent and TLI velocity differences are presented in Table 3.2. 
Close agreement was obtained between the results of the sensed and total 
velocity analyses, so that inclusion of the results for the total was 
unnecessary. Compensated sensed velocity differences (G&N - S-IVBIU, where 
the G&N data has been compensated with the derived IMU errors) are pre- 
sented in Figures 3-7 through 3-12. 

3.1.2 Ascent and TLI Error Determination 

The Apollo 13 G&C system accuracy analysis was based upon the deter- 
mination of a common set of errors which resulted in small residuals for 
both the boost to orbit phase and the translunar insertion phase. The 
analysis is accomplished with the aid of a Kalman Filter which solves 
for a "best" set of IMU errors for minimizing the velocity differences 
in a least squares sense. Several constraints were imposed on the 
errors used. The bias values for accelerometers (Table 3.3) and gyros 
were forcea to be in close agreement with inflight determined values 
and the other error terms were chosen to agree favorably with preflight 
calibration histories. Due to various physical factors such as actual 
parameter shifts during the boost phase and degradation of the reference 
data between the two flight phases (2.4 hours of drift between ascent 
and TLI) it was again recognized that all of the above conditions could 
not be met at all times. Based on engineering judgement, the approach 
pursued was to seek two sets of error sources with bounded variations 
lc )- The err <>r terms derived for the sensed analyses are presented 
in Table 3.2, and using these values, the G&N corrected trajectories 
fit the respective external measurement (S-IVBIU) trajectories. The maxi- 
mum deviation between the derived ascent and TLI error sources was O.860-. 
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The derived boost and TLI values for NBDX (X Gyro Constant Drift 
Rate) exceeded the lo (2 meru) instrument stability criteria. The de- 
rived values were 3.14 meru (ascent) and 2.35 meru (TLI). These 
represent 1.57a and 1.18a values respectively for boost and TLI. Pre- 
flight data obtained at Cape Kennedy, from 3 July 1969 to 3 April 1970, 
revealed a pronounced negative trend which peaked out at -3.3 meru on 
6 January 1970. Overall NBDX Cape test results are somewhat erratic. 
Following the referenced negative peak, NBDX results varied considerably. 
For the subsequent four calibrations, the quantity began trending positive 
and the last calibration value on 3 April 1970 was 0.5 meru. The CM com- 
puter erasable memory compensation value for NBDX was -0.7 meru. If the 
term continued to trend positively, it is probable that an effective error 
on the order of 2-3 meru did exist during ascent and Til. Consequently, 
it is understandable that the derived values for NBDX did exceed the la 
instrument stability criteria. The apparently large shift in drift prob- 
ably reflects a compensation error rather than Instrument degradation. 

All other error sources were within lo limits. 

3.1.3 Quick Look Evaluation 

It is worthwhile to point out that a technique has been developed 
for taking a gross "quick look" at IMU performance without recovering 
individual instrument errors. This was done for the Apollo 13 CM IMU 
by comparing actual ascent cut-off state errors with standard deviations 
of these errors. The standard deviations used for this comparison were 
obtained by integrating an ensemble of lo IMU instrument errors along 
the ascent trajectory. Resultant lo state errors formed trajectory 
bounds from which IMU performance could be gauged in the system sense; 
i.e., from which it could be determined whether or not overall IMU per- 
formance as a navigation Instrument was within 1 0 bounds. This gross 
"quick look" technique provided an additional confidence factor in the 
derived NBDX error source values discussed. Inasmuch as It demonstrated 
that IMU performance was not within specifications during ascent. A 
comparison of actual with lo velocity errors Is presented as follows 
(at t * 752 seconds GET): 
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X 

Y 

Z 

AX 

AY 

aZ 


(Ft)* 

(Ft) 

ifti 

(Ft/Sec) 

(Ft/Sec) 

(Ft/Sec) 

lo : 

1434 

27630 

3542 

4.46 

61.89 

9.88 

A Actual : 

3555.8 

32220.6 

1316 

-5.41 

75.51 

1.02 


*The position and velocity values are totals 
(u lits as noted) . 

3.1.4 Impact of Power Supply Degradation on CM PIPA Behavior 

Immediately following the SM LOX tank incident at 55:54:53 GET, 
the CM downlink recorded PIPA pulses which resembled the effect of venting. 
However, efforts to reconcile these data with doppler radar measurements 
met with failure. Consequently it was theorized that the actual pheno- 
menon being observed might be a bias shift in the PIPA's due to power 
supply transients, rather than accelerations due to venting. To support 
this hypothesis, an effort was made to correlate power supply transients 
with observed PIPA data. 

Figure 3-13 is a plot of the 28VDC Main A, 120VDC PIPA supply and 
the 3.2 KC 28V power supplies at the time of the LC^ tank event at 55:54:53 
GET. The 120V PIPA supply and the 3.2 KC 28V supply begar. to degrade when 
the 28V0C Main A supply reached 26.3 volts. Figure 3.14 is a plot of the 
accumulated thrust velocity (V^, V^y, V^) indicated by the PIPA's through 
the time of the voltage transient; superimposed on this is a plot of MSFN 
doppler radar residuals for the same time period. The change in PIPA 
velocities can be directly correlated to the degraded voltage. The reason 
for the different PIPA responses (i.e., XPIPA does not show response to the 
sharp voltage transient at 56:00 and 56:03 as do Y and Z) is unknown. 
Referring to the plot of the Bermuda/Goldstone tracking station doppler 
residuals through the period of Interest, no appreciable change is noted 
between 55:57:56 and 56:03:10, the region of the 120V PIPA supply transient. 
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Shifts in the doppler residuals near 55:55:00 and 56:03:30 resemble velo- 
city changes indicated by the X PIPA which is pointed principally along 
the LOS to Earth. The trend shown by the X PIPA however, is not reflected 
in the doppler data, thus indicating that X PIPA velocity indication- were 
erroneous. 

Following the event, the CSM IMU was powered down (including instrument 

heaters) at 58:40 GET and remained down until 140:10 GET. At IMU power up, 

2 

a shift of 1.64 cm/sec was noted on the Z PIPA and as a result a bias update 
was performed prior to entry, lable 3.3 is a summary of the PIPA bias 
before and after the event. 

3.2 CM OPTICAL NAVIGATION 

Prior to the SM LOX tank incident, two batches of navigation sight- 
ings were conducted. The data obtained from these sightings is listed 
in Table 3.4. These data were analyzed for the purpose of estimating 
sextant trunnion bias and earth horizon bias. 

Two basic options were available for analysis of these data. In 
Option 1, only trunnion angle data Is processed. This option assumes 
that the astronaut was successful in locating the substellar "point" - 
l.e., In placing his horizon llne-of-slght on the line connecting the 
star with the earth's center. Option two utilizes glmbal and shaft 
angle data for the purpose of correcting erroneous determinations of 
the substellar point. In addition to these options, batch one was pro- 
cessed using: 

Method 1: Horizon altitude only estimation, and 

Method 2: Horizon altitude and trunnion bias estimation, 
with no a priori estimate of these parameters 
(l.e., with Initial estimates weighted such that 
they had no effect on the solution). 
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The combination of method? one and two with options one and two 
yielded four cases. The output of these cases were used as the input 
(initialization values) for the processing of batch two. Table 3.5 
presents these combined results. Option 1 results from batch 1 were 
used only with option 1 runs for batch 2. Likewise, option 2 results 
from batch 1 were used only with option 2 runs for batch 2. 

The results indicate that the actual horizon altitude was about 18 
kilometers. However, due to small sighting errors the astronaut was 
actually sighting on an effective altitude of about 10 kilometers. 
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Table 3.3 CM IMU PIPA BIASES 


GET 

2 

X Bias (cn/sec ) 

2 

Y Bias (cm/sec ) 

Z Bias (on/sec 1 

10:10 

-0.20 

-0.21 

0.0 

12:25 

-0.20 

-0.20 

0.0 

13:05 

-0.19 

-0.21 

0.0 

14:35 

-0.20 

-0.20 

0.0 

16:35 

-0.20 

-0.20 

0.0 

35:09 

-0.20 

-0.21 

0.01 

42:48 

-0.20 

-0.20 

0.0 

57:16 

-0.21 

-0.12 

-0.01 

141 :00 

-0.18 

-0.16 

-1.66 
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Table 3.5 COMBINED RESULTS OF PROCESSING 
BATCHES 1 AND 2 OF P23 DATA 


Run 1 
Method 2 

b 

°h 

°b 

p 


Run 2 
Method 2 


Run 3 
Method 1 

b 


°h 


°b 

p 


Run 4 
Method 1 


b 


°h 


°b 

p 


Initialization 

Option 1 Results 

Option 2 Results 

10.61 km 

10.13 


0.00952 mrad 

-0.02236 


4.88 km 

4.57 


0.07975 mrad 

0.05386 


0.27406 

0.12275 



17.42 

17.55 

-0.02827 

-0.03093 

4.88 

4.57 

0.07975 

0.05386 

0.27406 

0.12276 


9.61 

9.86 

-0.05236 

-0.05236 

4.69 

4.53 

IQ' 5 

10" 5 

0.0 

-0.00001 


17.04 

17.35 

-0.05236 

-0.05236 

4.69 

4.53 

IQ' 5 

10" 5 

0.0 

-0.00001 


where , 

h = estimate of altitude bias 
b = estimate of trunnion angle bias 

o h = standard deviation of estimate of altitude bias 
from covariance matrix 

° b ■ standard deviation of estimate of trunnion bias 
from covariance matrix 

p ■ correlation coefficient of estimates of altitude and 
trunnion bias from covariance matrix 
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4.0 LM SYSTEMS 


4.1 LM INERTIAL MEASUREMENT UNIT 

4.1.1 Trar.searth Injection Burn: Cut-Off Velocity Errors 

4. 1.1.1 Velocity Errors 

Velocity errors generated by the IMU during the TEI burn were in- 
vestigated. To do this, "best" ignition and cut-off vectors were deter- 
mined from free flight data. Three such vectors were determined for the 
time of ignition; one was determined for the time of cut-off. These 
vectors are presented in Table 4.1. The ignition vectors were extra- 
polated to the time of cut-off using accelerometer data from the LGC 
downlink. These extrapolated vectors may be thought of as LGC "best 
guesses" as to the vehicle state at time of cut-off. They were differ- 
enced with the free flight cut-off vector (assumed correct) to obtain 
estimates of the IMU induced trajectory errors. Only the velocity 
components were calculated, and the results of these calculations are 
presented below for each ignition vector and each of three coordinate 
frames. The coordinate frames chosen were: BRCS (Basic Reference 

Coordinates); IMU, and; local. All reference frames were moon centered. 

The local frame is defined to agree with MIT's convention, as follows: 

X: completes the right hand set 

Y: along the negative angular momentum vector {-RxV_) 

Z: along the negative of the radius vector (-R) 

A V aV aV 

("Best" - NBI.X343) (GYMX289 - N BEX343) (HAWX300 - NBEX343) 

XY_ZXY_ZXY£ 

(+3.85; +3.02; +2. i2) (+2.92; -1.21; +9.62) (+3.53; -0.14; +8.44) 

(-3.77; -3.83; -1 .07)(-2.72; -2.92; -9.31) (-3.23; -3.53; -7.80) 

(+5.46; +0.41; >0.56) (+5.10; +8.72; -0.72) (+5.70; +7.12; -0.66) 

5.48ft/sec 10.13 ft/sec 9.15 ft/sec 


BRCS 
IMU 
LOCAL 
Magni tude 









* 


The GYMX 289 and HAWX 300 ignition vectors yield similar velocity 
errors whereas errors for the "best" ignition vector differ from these 
markedly. The first two vectors were obtained by the RTCC in real time. 
GYMX 289 was the vector with which the TEI burn was targeted. However, 
both vectors contain only data obtained prior to perilune. The "best" 
vector was obtained at TRW by fitting data on both sides of perilune. 

It provided much more satisfactory range and doppler residuals after 
perilune and consequently has been chosen as the preferred vector for 
purposes of this report. The disagreement between the preferred vector 
a.id the other two vectors lie almost entirely in the out of plane 
("Y", in local coordinates) direction; this is the direction in which 
velocity is least accurately determined by ground based radar (and, 
therefore, the one in which the greatest uncertainties are experienced). 

4 . 1 . 1 . 2 F .rror Sources 

The principal sources of velocity error are likely to have been 
IMU misalignments (o.7 deqrees/axis , due to l<r drifts alone) and dis- 
crepancies in the ignition and cut-off vectors obtained from free flight 
data. These free flight vector uncertainties are not known explicitly. 
However, the ignition vector is credited with very good accuracy; on the 
order of .1 ft/sec, out-of-plane, and less than .1 in the other two 
axes. The cut-off vector is suspected of greater errors, but these may 
well be less than one ft/sec (total). With that hypothesis, it was 
concluded that the majority of the error could be attributed to IMU 
misalignments. Since alignment errors are perpendicular to the velo- 
city gained, it was of interest to determine how much of the observed 
error satisfied that condition. This was done as follows: 

V p (parallel to VG) = (V) • (unit VG) * -1.47 ft/sec 

Vp = V p (unit VG) = [-1.27; +0.73; -0.16), In IMU 
coordinates 

(normal to VG) * V - V^ • [-2.50; -4.56; -0.91], 

IMU coordinates 

V n * 5.27 ft/sec 
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Most of the error was perpendicular to VG, tending to support the 
argument for IMU misalignments as the principal factor. The -1.45 ft/sec 
of error parallel to VG cannot be attributed to IMU misalignments and since 

other IMU contributions can reasonably be expected to have been much less 
than this, it seems probable that the free flight cut-off vector contained 
a -1.47 ft/sec error along VG (and an unknown, but presumaDly small com- 
ponent perpendicular to it). 

4.1.2 IMU Misalignments at Time of TEI Burn 

Three approaches were taken to determine the IMU alignment errors 
which existed at the time of the burn. These are described in the 
following paragraphs. 

4. 1.2.1 Sun Check 

The crew sun check (discussed in Section 4.4.1) performed at 73:47 
GET provided information concerning IMU misalignments perpendicular to 
the line-of-sight. Unfortunately, there were Insufficient data to deter- 
mine the direction of those errors which were observed. However, it 
was determined that the magnitude of observable errors was approximately 
0.51 degrees. This provided a useful method of testing the validity 
of other estimates. 

4. 1.2. 2 Drift Check 

An effort was made to determine alignment errors from gyro drift 
data. These data included LM prelaunch and CM inflight measurements. 
Basically, it was assumed that LM IMU error at any given time was the sum 
of: 

o CM IMU drift prior to the docked alignment, plus 

o LM IMU drift subsequent to the docked alignment. 

The uncertainty in the estimate thus obtained was taken to be the RSS 
of CM and LM initialization errors. Uncertainties in the drift rates 
were Ignored. The point of this was to determine whether misalignments 
could have been reliably deduced from existing Instrument performance 
data and eliminated prior to the midcourse burn. Table 4.2 presents 
the results of these estimates. 
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From the table it can be determined that the magnitude of misalign- 
ment estimated from drift data is 0.60 degrees, with a la uncertainty 
due to initialization of +.28. Almost all of this would have been 
visible to the sun check. Thus, the magnitude of the predicted error 
agrees rather well with that observed in the sun check, and the dis- 
agreement lies well within the bounds of predicted uncertainty. 

4. 1.2. 3 Velocity Error Check 

Once the velocity error at TEI cut-off was determined, it was 
possible to determine in two steps a unique set of alignment errors which 
satisfied the constraint that the sum of the squares of the error angles 
was minimized. Step one consisted of determining that component of 
velocity error perpendicular to the velocity gained vector. This was 
done because IMU misalignments can be shown to produce only perpen- 
dicular errors. The procedure and results were described in Section 
4. 1.1. 2. The second step consisted of solving for a "minimum" set of 
error angles. The constraint equations are: 

Minimize * i = ! t ! 2 , subject to 

(VG) x (i) - AV , where 

*~~TI 

j_: the vector of error angles 

vf: the velocity gained in the burn 

j,V : that component of burn velocity error perpendicular 

n to VG. 

The solution obtained from these constraints is: 

i-(vi) J UV «») 

Applying this relationship to the quantities generated in Section 4.1.1 
yields: 


= [-.062; -.035; +.344] degrees 
’;!= .351 degrees 
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It must be understood that, assuming a correct measure of aV^, these are 
the errors which generated MT (moreover, most of this would have been 
visible in the sun check). However, these errors do not necessarily 
reflect total IMU misalignment. Platform misalignments about VG would 
generate no velocity error. Consequently, that component of total IMU 
misalignment which was about. VG cannot be determined by an analysis of 
aV^. It is to be expected that total IMU misalignment magnitude would 
be no less than .351 and would probably be greater. The discrepancy 
would arise from the irresolvable component about VG. 

A pattern of reasonably good agreement emerged from the above error 
estimates as evidenced in the summary below. 

Source of Estimate <t> . a i a t 

_i _x _y. _z 

Sun Check 

Drift Check -0.26 +0.24 -0.07 +0.14 +0.54 +0.04 

Velocity Check -0.06 -0.04 +0.34 


0.51 

0.60 +0.28 
0.35 


Agreement between the three methods in the Z component, and in magni- 
tude, is fairly good. Also, agreement in magnitude and in the "X" and "Y" 
estimates lies within the 1c uncertainty band. 


4.1.3 PIPA Bias and IRIG Drift 

LM PIPA biases were stable throughout the periods the LM IMU was 
activated. The LM IMU was turned on for the first time shortly after 
the CSM incident and remained powered up until after the TEI burn. Heater 
power was always maintained in the system. The system was activated 
for the second time at 134 hours GET in preparation for the MCC 4 maneu- 
ver. Samples of PIPA bias during the periods the system was activated 
Indicated the following mean values and data variations. 
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PIPA 

Number of 
Samples 

Mean Value 

Sample* 

Standard Deviation 

X 

19 

1.494 cm/sec 2 

0.014 cm/sec 2 

Y 

19 

-1.427 cm/sec 2 

0.038 cm/sec 2 

Z 

19 

1.573 cm/sec 2 

0.011 cm/sec 2 


♦Design Specification Uncertainty (1 j) = 0.2 cm/sec 2 . 


Insufficient data were available to evaluate IMU IRIG drift on this 
mission. However, based on the small values of LM IMU misalignment at 
the time of the TEI bum determined from the error separation study in 
Section 4. 1.2. 3, it can be deduced that the gyro drifts were small and 
easily within the 3c design uncertainty of 0.09 deg/hr. 

4.2 LM DIGITAL AUTOPILOT 

An analysis of the DAP orirol functions during the Apollo 13 mission 
was performed to verify proper DAP performance. The following items were 
considered: 

o The MCC 2 DPS free return maneuver 
o The TEI (peri cynthion + 2 hour) DPS maneuver 

o The maneuver to the LM PGNCS PTC attitude. 

o The attitude hold capability of the LM DAP in the 
LM/CSM and LM/CM configuration. 

o The automatic maneuver capability of the LM DAP in 
the LM/CSM configuration. 

o The manual maneuver capability of the LM DAP in the 
LM/CSM and LM/CM configuration. 

4.2.1 MCC 2 DPS Free Return Maneuver Analysis 

A four jet, 8.01 second ullage was initiated prior to the MCC 2 - 
DPS burn at 61:29:36.06 GET. DPS ignition occurred at 61:29:42.84 GET. 
Prior to the DPS burn the U-V rotational jets (i.e., X translation jets) 
were not manually inhibited as is usually the procedure. Manual throttle- 
up from 11.9* to 36.8* occurred at approximately 61:29:50 GET. Program 
sequencing prior to and during the bum was nominal. Manual throttle-up 
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to 36.8% was nominal and required aV of approximately 38 ft/sec was 
achieved. The burn residuals were 0.2, 0, and 0.3 ft/sec for the X, 

Y, and Z components, respectively. No manual nulling of these residuals 
was necessary. Table 4.3 shows the magnitudes of the maximum estimated 
rates and the rate gyro signal maximum values during this burn. The 
low magnitudes indicate good burn performance and no discernible 
sloshing effects. Little or no slosh would be expected since the APS 
and SPS were fully loaded during the burn and the DPS was fully loaded 
at the start of the bum. Table 4.4 shows the magnitudes of the maximum 
attitude errors and rate errors obtained during this burn. Figures 4-1, 
4-2, and 4-3 show the phase plane plots for the P, U, and V axes, re- 
spectively, for the MCC 2 burn (Relationship between P, U, V and X, Y, Z 
axes is shown in nomenclature). These plots indicate nominal DAP perfor- 
mance. Table 4.5 shows the RCS fuel consumption required to maintain 
attitude control during this burn. The total RCS propellant (excluding 
ullage) required during MCC 2 was 12.90 lbs. Figures 4-4 and 4-5 show the 
pitch and roll GDA positions throughout this burn and indicate satisfactory 
GDA performance. During the manual throttle-up from 11.9% to 36.8%, the 
effects of a transient due to the compliance of the DPS gimbal system can 
be seen. This compliance effect results In the GDA's driving in + pitch 
and a - roll direction. After relieving this compliance, the GDA movement 
Is as expected in trying to track the c.g. 

4.2.2 TEI DPS Maneuver Analysis 

Prior to the TE! hum, jets 6 and 14 were used for a two jet, 11.1 
second ullage which began at 79:27:23.43 GET. DPS ignition occurred 
at 79:27:38.30 GET. Prior to the burn, the U-V rotational jets were 
manually Inhibited as required by the preburn checxllst. A nominal 
manual throttle-up from 11.4% to 37.7% occurred at approximately 
79:27:45 GET. The automatic throttle-up to 93.9% occurred at approxi- 
mately 79:28:05 GET. The burn residuals were 1.0, 0.3, and 0 ft/sec for 
the X, Y, and Z components, respectively. No manual nulling of these 
residuals was performed. Table 4.6 shows the magnitudes of the maximum 
estimated rates and the rate gyro signals maximum values during this burn. 
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The values indicate nominal performance and very smull slosh effects, 
as would be expi .ed with the propellant loadings that were present, 
i.e., APS full and DPS and SPS nearly full. Table 4.7 presents the 
magnitude of the maximum attitude errors and rate error which occurred 
during this burn. Figures 4-6, 4-7, and 4-8 show the phase plane plots 
for the P, U, and V axes, respectively. The time period plotted begins 
10 seconds before ullage and extends 20 to 30 seconds beyond the start 
of automatic throttle-up. The plots indicate nominal DAP performance. 

The U and V axes plots indicate two large excursions beyond the dead- 
band. In each case, the first excursion is smaller and is associated 
with manual throttle-up. The second excursion is associated with auto- 
matic throttle-up. Table 4.5 presents the RCS fuel consumption required 
to maintain attitude control during the burn. The total RCS propellant, 
(excluding ullage) was 4.90 lbs (all expended in P-axis control). 

Figures 4-9 and 4-10 show the pitch and roll GDA positions from DPS 
ignition through the automatic throttle-up maneuver and indicate nominal 
GDA performance. Some concern was expressed about the fact that the roll 
GOA drove approximately -1.3 degrees from its initial oosition at the 
start of the TEI burn. This appears to be nominal behavior caused by 
relieving compliance and perhaps correcting for some mistrim, The 
following facts are known. 

a) The preferred direction for compliance Is in a "+ pitch" 
and a "- roll" direction. 

b) Compliance appears to be highly non-linear and transient 
effects cannot always be observed. 

c) Previous simulations using the now- inoperative MSC bit- 
by-bit simulation have shown similar transient compli- 
ance effects on all tests. 

d) For the MCC 2 bum the pitch GDA drove +1.05 degrees 
during throttle-up from 11. 9% to 36.8% and the roll 
GDA drove -0.75 degree during throttle-up from 11.9% 
to 36.8%. These effects were apparently due to the 
effects of compliance. 

e) At the end of the MCC 2 bum both the pitch and roll 
GDA's were driving in a positive direction In trying 
to track the c.g. 
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f) For the TEI burn the pitch 6DA drove +0.435 degree 
during throttle-up from 11.4% to 37.7% and the roll 
GDA drove -1.335 degrees during this throttle-up 
period. 

g) The root- sum- square GOA change due to compliance was 
approximately 1.29 degrees for the MCC 2 burn and 
1.402 degrees for the TEI burn. 

h) For the TEI burn the pitch GDA showed no noticeable 
change during throttle-up from 37.7% to 93.9%. The 
roll GDA drove -1.245 degrees during this throttle-up 
period. The root-sum-square GDA change of 1.245 degrees 
agrees closely with the previously discussed values of 

1 . 29 degrees and 1 . 402 degrees . 

Evaluation of the facts discussed above leads to the conclusion that the 

GDA behavior for both the MCC 2 burn and the TEI burn was nominal. 

4.2.3 Maneuver To PTC Attitude 


The difficulties encountered in maneuvering to the PTC attitude 
subsequent to the TEI burn prompted an investigation of the torques 
associated with each LM reaction control jet and with the various 
rotational and translational firing policies available in the LGC. 
This investigation included the effects of Impingement forces on the 
plume deflectors. 

At the time of the maneuver, the best estimate of LM and CSM 
weights were 25666.2 lbs and 62489.7 lbs, respectively, or a total 
weight of 88155.9 lbs. This value was referred to in the SODB 
(Reference 2) which defined the combined vehicle c.g. location in the 
CSM body frame as follows: « +1049.38 Inches; * +2.49 Inches; 

and « +3.59 Inches. These values were obtained by linear inter- 
polation of the values recorded In the SODB. These data were used to 
transfcnn the c.g. location into LM body axes. 


4-9 



The SODB also defined the locations of the points at which the re- 
action control jet forces are applies. The location of the impingement 
forces on the plume deflectors was defined by Reference 3. Based on the 
location of each of the sixteen reaction jets with respect to the c.g. and 
the thrust vector, the torque for each jet was calculated. Using the data 
from Reference 3, torques were calculated for the impingement forces result- 
ing fro,. i each of the downward firing thrusters. 

The GSOP (Reference 4) gives the iet policies to be used in the 
various translation and U-V rotational maneuvers that can be performed. 
Tables 4.8 and 4.9 indicate the torques which result from the use of 
each of these policies. In addition, there are possible jet policies 
associated with Q and R axes rotations in the minimum impulse mode. 

These are not defined in the GSOP. However, the various i.olicies have 
been determined and the torques resulting from each policy are defined 
in Table 4.10. 

The maneuvering problem reported by the astronauts primarily was 
a cross-coupling effect between the pitch (Y or Q) axis and the roll 
(Z or R) axis. No downlink data is available for the period of the 
maneuver. However, certain assumptions can be made. 

a) All jets were functioning; this implies that primary 
jet-policies were in effect at all times. 

b) The maneuver was conducted manually using either: 

1) Y-axis and Z-axis translation commands through 
the TTCA; or 

2) Three-axis (X, Y, and Z) rotational commands 
through the ACA. 

Examination of Table 4.8 reveals chat Y-axis and Z-axis translation com- 
mands produce only minor cross-coupling around the yaw axis. Pitch-roll 
cross-coupling is non-existent. Examination of X-axis rotational policies 
in Table 4.9 reveals that primary jet-policies, using either four jets 
or two, produce only X-axis rotations. Data for Y-axis and Z-axis rota- 
tions are shown In Table 4.10. The use of four jets for any rotational 
maneuver produces very minor cross-coupling. However, the GSOP states 
that pitch and roll maneuvers In the minimum Impulse mode will be 



accomplished with two jets. Examination of the jet-policies which use 
two jets reveals cross-coupling torques that have the same magnitude as 
the primary torques. This occurs because each of these jet-policies 
employs one downward firing jc'. The resulting impingement force acts 
equally about the Y-axis and the Z-axis using the long X-axis moment 
arm. The result is to reduce the primary torque by a factor of two 
and produce the large cross-coupling force. 

Although data for this time period is lacking, examination of data 
for a later period revealed that these rotational firings took place. 

Data indicated a +R (+Z) rotational command at 137:22.49.516 and -0 (- v ) 
rotational command at 137:23:01.616. Examination of jet firings data 
showed corresponding 60 millisecond firings of jets 5 and 10 at 137:22:50.434 
and of jets 10 and 13 at 137:23:01.234. Analysis has shown that jets 5 
and 10 produce a primary +R torque and jets 10 and 13 produce a primary -Q 
torque. Since two jet Y and Z axes torques occurred at these times, it 
is not unreasonable to infer that they also occurred during the maneuver 
to PTC attitude and that they produced not only the primary torque but a 
cross-coupling torque of equal magnitude as well. If the THC was used 
to perform maneuvers in the LM/CSM docked configuration, no unexpected 
problems should be encountered if the crew is familiar with this type 
of operation. However, use of the RHC or the MIC can produce results 
which would be very unexpected. 

It is recognized that the above analysis represents only one of a number 
of possible explanations of the control difficulties encountered. Another 
probable source of difficulty is the fact that the FDAI was powered down 
during the maneuver. Without that source of attitude information, it was 
necessary to monitor attitude by observing glmbal angles displayed on the 
DSKY. Because the spacecraft yaw axis was not coincident with that of 
the platform, either a pitch or a roll command would cause a change in both 
of the corresponding gimbal angle displays. That this was a factor is sup- 
ported by the fact that, after the vehicle's and platform's corresponding 
axes were brought into closer alignment, passive thermal control was estab- 
lished satisfactorily. 
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4.2.4 DAP Control of the LM/CSM Configuration 


DAP attitude control was evaluated for periods of attitude hold, 
and for periods of automatic and manual maneuvers. The available data 
were for a 5 degree deadband. Figures 4-11, 4-12 and 4-13 (P, U, and 
V axes, respectively) present both periods of auto maneuver and attitude 
hold. The auto maneuver period was 61:25:17 to 61:25:20 GET. The atti- 
tude hold period began at 61:25:21 GET and was plotted for approximately 
2 minutes. This period precedes the MCC 2 maneuver and nominal DAP 
performance is indicated. An auto maneuver prior to the TEI burn is 
presented in Figures 4-14, 4-15, 4-16 (P, U, and V axes, respectively). 
The maneuver began at 79:21:43 and was terminated 3 seconds later. 

Nominal DAP performance is indicated. A five minute period of attitude 
hold subsequent to this maneuver is plotted in Figures 4-17, 4-18, and 
4-19 (P, U, and V axes, respectively). The period covered is 79:22:00 
to 79:27:00. The U and V axes phase Dianes are plotted for the entire 
time period. The P axis was functioning in a much faster duty cycle of 
which one complete cycle has been plotted. Manual maneuvers were per- 
formed in the minimum impulse mode where the desired CDU angles are 
set equal to the actual CDU angles. Since no attitude errors exist, 
it is not possible to generate a meaningful phase plane plot. However, 
it was determined that the rate errors did not exceed the deadband 
limits. 

4.2.5 DAP Control of the LM/CM Configuration 

DAP attitude control was evaluated for a time period that contained 
both attitude hold and a manual maneuver. The P, U, and V axes phase 
plane plots are presented in Figures 4-20, 4-21 and 4-22, respectively. 
The attitude hold period presented began at 140:49:12 GET and terminated 
at 140:51:25 GET at which time the manual maneuver began. Since the 
manual maneuver was performed in the minimum impulse mode, the attitude 
error remained approximately zero as shown by the three figures. 
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4.3 LM ABORT GUIDANCE SYSTEM 






Investigation of AGS performance was conducted with the objectives 


o Determining the stability of AGS sensor static errors 
(accelerometer bias; gyro static, or bias drift). 

o Determining AGS sensor dynamic errors from comparisons 
of AGS and PGNCS measurements during the TEI burn. 

o Comparing AGS and PGNCS measurements of velocity gained 
and vehicle attitude during midcourse burns. 


After the Service Module LOX tank incident, four burns were made to 
return the spacecraft to earth as quickly and safely as possible. The 
first and second burns were made under PGNCS control with the LM DPS; 
the AGS was powered up for the second burn and used in the "back-up" 
mode in case of PGNCS failure. The AGS was used to perform the last 
two burns, a DPS burn of 7.8 ft/sec and an RCS burn of 2.8 ft/sec. 

4.3.1 Burn Analysis 

TEI (AGS in Follow-Up) 

The second LM DPS burn was initiated at approximately 79:27:41 
following an X-axis RCS burn for ullage. This burn was performed with 
PGNCS control and the AGS in follow-up. Although the AGS was not tar- 
geted for the maneuver, its indication of sensed velocity and attitude 
were available to confirm success of the burn. After AGS power-up, 
a Body Axis Align was performed. Throughout most of the burn, the astronauts 
monitored "X" sensed velocity in body coordinates on the DEDA and PGNCS sensed 
velocity on the DSKY. These velocities are plotted in Figure 4-23 and 
indicate good agreement between the AGS and PGNCS. 

Following the burn, AGS power-down (including heaters off) began at 
approximatley 79:51:00. 
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MCC 3 (AGS Controlled - DPS) 

This burn, a small maneuver to increase the entry flight path 
angle, was performed under AGS control. Plots of sensed velocities 
are shown in Figure 4-24. Since the burn was insensitive to burn time 
cut-off errors and attitude errors, the primary burn rule was to avoid 
excessive rates (more than 10 leg/sec) about any axis. 

System power-up for the burn was initiated at approximately 
104:40:00. At approximately 104:58:20, the ASA temperature had risen 
above 115 degrees F and the AGS was activated. 

DPS burn time was voiced up as 15 sec although the astronauts were 
told to shut-down after 14 sec to avoid an overburn, which would require 
an RCS trim, impinging on the command module. 

Figures 4-25, 4-26, and 4-27 show the attitude rates and attitude 
errors. Rate and attitude error data indicate no appreciable disturb- 
ances since all three channels appeared to be in normal limit cycle 
operation. 

Because of the planned early shutdown, the maneuver was a slight 
underburn. An RCS trim burn was performed approximately one minute 
later, increasing the burn aV by 0.2 ft/sec to 7.8 ft/sec. Subsequently, 
the AGS was powered-down including heater power off. 

MCC 4 (AGS Controlled - RCS) 

PGNCS and AGS were powered-up at approximately 136 hours and the 
maneuver to burn attitude for this midcourse correction was with the 
PGNCS . 

In Figures 4-29, 4-30, and 4-31 comparisons between the AGS and 
PGNCS indicated attitudes are made for the period of manual maneuver- 
ing to acquire the appropriate spacecraft attitude for MCC 4. (An 
AGS to PGNCS align preceded the re-orientation.) Following acquisition of 
the burn attitude, Guidance/Control was switched to AGS. 
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The burn was initiated at approximately 137:39:52 with the X-axis RCS. 
The DSKY was used «.o observe the accumulated velocity. After about 23 
seconds, the RCS was turned off after achieving a aV of 2.7 ft/sec. A 
partial trim was effected 15 sec later, bringing the total aV to 2.8 ft/sec. 
This was short of the desired value of 3.1 ft/sec, but was indicated as 
acceptable by the ground. The aV's are shown in Figure 4-28. 

No attitude control problems were noted during the burn, with manual 

pitch and roll provided by the TTCA and AGS yaw control by the AGS at:itude 
hold mode as seen in Figure 4-32. 

4.3.2 Sensor Performance (ASA 023) 

4. 3. 2.1 Accelerometer Errors 

Free Flight Performance - Free flight accelerometer static bias data 
were obtained for two time periods, pre-TEI and post-MCC 4. Bias esti- 
mates obtained from these data appear in Table 4.11. Table 4.12 compares 
the shifts undergone by this parameter after prelaunch calibration with 
standard deviations of this shift. The AGS Capability Estimate one sigma 
values result from data samples taken from a number of AGS systems; as such, 
they are representative of general AGS performance. The table reveals that 
static bias performance was easily within 1? of AGS Capability Estimate. 

This is excellent performance in view of the fact that the ASA dropped 

o 

23 degrees F below the specified minimum (60 F) during the 24 hour period 
in which it was shutdown prior to the TEI burn. 

Powered Flight Performance - The TEI burn was the only one performed 
during Apollo Mission 13 of sufficient duration and thrust to permit 
estimation of AGS dynamic sensor performance. Figures 4-33, 4-34, and 4-35 
show the accumulated sensed velocity along the body axes. The aV magnitude 
was about 860 ft/sec. 

Accelerometer errors were found from AGS/PGNCS sensed velocity 
comparisons. The PGNCS velocities were compensated for known static 


4-15 





bias errors, interpolated to AGS times, transformed to body coordinates, 
and subtracted from the corresponding AGS velocities. The differences 
are shown i.. Figures 4-36, 4-37, and 4-38. Since PGNCS gimbal angles 
were used to transform the velocities, these differences reflect only 
accelerometer errors (plus noise due to AGS downlink and PGNCS gimbal angle 
quantizations) . 

The estimates of dynamic and total accelerometer error are listed in 
Table 4.13. The total errors were derived with a digital computer pro- 
gram which determines a set of "best" AGS errors for the purpose of mini- 
mizing the AGS-PGNCS velocity residuals in a least squares sense. The 
dynamic errors were obtained by subtracting the static bias (measured over 
a 15 minute interval just before the burn) from the totals. Table 4.15 
presents ratios of the shifts (from prelaunch values) in these parameters 
to shi-f + s predicted by the ASA 023 error model. The table shows that all 
of these errors were within the error model 2o limits. Plots of the velocity 
residuals, compensated for the recovered errors, are shown in Figures 4-39, 
4-40, and 4-41. 

4. 3. 2. 2 Gyro Errors 

Free Flight Performance - Free flight gyro static drift data were 
obtained during two time periods, pre-TEI and post-MCC 4 (the interval ' *om 
140:29 to 140:51 GET). Drift estimates obtained from these data appear n 
Table 4.11. X channel data were noisy and interpretation was correspondingly 
difficult. Thus it was necessary to specify a range of possible performance 
values. Instrument performance was within lo AGS Capability Estimate 
tolerances, so that it may be said that the static drift performance was 
quite satisfactory. 

Powered Flight Performance - Gyro error was obtained from AGS/PGNCS 
attitude differences obtained just before and during the bum. These 
differences are shown in Figure 4-42, 4-43, and 4-44. Because of 
relatively constant inertial attitude and thrust level, the effects 
of dynamic drift, scale factor ernor, input axis misalignment, and mass 
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unbalance cannot be separated. However, since the burn was performed 
in the "attitude hold" mode, the effect of a significant scale factor 
error or input axis misalignment would have been negligible. Because 
of the low thrust level (about 0.1 g), mass unbalance probably con- 
tributed very little error. Thus most of the error was attributed to 
gyro static and dynamic drift. The total errors were determined by 
measuring the slopes of the attitude differences during the burn. 

The dynamic error was obtained by subtracting the static values from 
the totals. The dynamic and total errors are presented in Table 4.14. 

Listed in Table 4.15 are ratios of the shifts in measured parameters 
to their la values. None of the errors exceeded 2o with respect to ASA 023 
prelaunch performance. 

4.4 LM OPTICAL ALIGNMENT CHECKS 

4.4.1 Sun Check for TEI Alignment 

Shortly after the service modul° LOX tank incident, the crew per- 
formed a docked LM 1MU alignment to the CM IMU. At 73:47 GET a check 
of this alignment was performed. At the time of the check, visibility 
through the optical instruments was extremely poor due to an agregation 
of debris surrounding the spacecraft. In addition it was desired that 
maneuvering be held to a minimum to conserve Ru> propellant. Conse- 
quently the decision was made to limit the initial check to one celestial 
body. The sun was chosen because of the visibil'-y conditions. The 
check consisted of pointing the AOT line-of-sight toward the LGC calculated 
sun direction and then noting the approximate magnitude of the offset. 

The LGC was placed in P52 and routine R52 (Auto Optics Positioning) was 
used to establish the LGC-estimated pointing vector. The FDAI attitude 
error displays were used to hold this spacecraft attitude while the 
position of the sun in the field of view was determined. It was deter- 
mined on the ground that the solar disc would subtend approximately 
one-half of a degree In the AOT field of view. From that fact and the 
observed sun Image the crew determined an Image (and therefore an IMU) 
misalignment of approximately one-half of a degree in an undetermined 
direction. Since this single target optical check could not resolve 
errors about the line-of-sight, no estimate of possible misalignments in 
that direction was obtained. The mission rule established for this 
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check was that any observed terror of less than one degree was tolerable. 

No effort was made to correct the error observed, or to determine the 
error about the line-of-slght, prior to the TEI burn. 

In conducting a single target alignment check such as this it is import- 
ant to select a line-of-sight which reveals those alignment errors of con- 
sequence in subsequent powered flight phases. More explicitly, it 
should be chosen to measure all misalignments having a significant impact 
on trajectory parameters of importance to a safe re-entry. One very 
crucial parameter is entry flight path angle, and during the deep space 
phases of transearth flight this is governed largely by inplane velocity 
perpendicular to the earth centered radius vector. The sun was an ex- 
cellent choice of targets for detecting misalignments as they effected 
Apollo 1 3 ' s post-TEI entry corridor. This is demonstrated by the 
results presented below: 

Let . = entry flight path angle 

V = velocity error in the direction which effects flight path ang 

i 

V = velocity gained in the TEI burn 

Unit (L0S)= the unit line-of-sight vector to the sun at the 
time of the sun check. 

*L0S s I,1U misa ^9 nment about the line-of-sight vector 

Hv “| 

Then ^ s #1 L Y- • (V X Unit LOSj J * -.088 

This quantity is unltless and serves to show that a one degree misalign- 
ment (for example) of the IMU about the line-of-sight vector would have 
produced only -.088 degrees of error In the entry angle of the post-TEI 
trajectory. Restated, entry flight path angle errors were extremely 
Insensitive to misalignments which the sun check could not resolve. 

Conversely, a was quite sensitive to those errors which could be 
observed by the sun check, as Illustrated below: 

It- * -10., where Is misalignment observed In 
n 

the sun check. 



Based on the above study, it is concluded that: 

o r he choice of the sun as the target for this optical check 
was an excellent one, and; 

o Compensative torquing of the IMU to eliminate the error ob- 
served in the sun check would have resulted in a much safer 
post-TEI trajectory. 

4.4.2 Sun/Moon Alignment Star Angle Difference 

During the time interval from 134:45 through 135:02 of ground elapsed 
time the Apollo 13 crew performed a P52 alignment of the LM IMU. The 
sun and moon (centers) were the optical targets. A star angle difference 
of -1.12 degrees was calculated by routine R54, indicating a very large 
astronaut sighting error. However, a postflight investigation shows 
the actual star angle difference to have been approximately 0.08 degree. 

The LGC has an algorithm and associateo ephemeris data stored within 
it for computing pointing vectors to various celestial bodies, Including 
the sun and moon. For the star angle difference calculations, this al- 
gorithm is used and the scalar product of the resulting vectors is 
computed. The arc cosine of this quantity forms the LGC's best estimate 
of the angle between the sighting vectors. The LGC also computes the 
angle between the measured vectors, and these two angles are differenced 
to obtain an error (stored - measured) for presentation to the crew. 

After the above mentioned Apollo 13 sightings, the crew display showed 
an angular error of -1.12 degrees. In order to check this value, the 
two angles and their differences were recomputed independently as 
described below. 

Procedural ly, the crew can make as many as five measurements of each 
vec^r. The IGC averages these measurements and transmits the two aver- 
age vectors on the downlink. These veciors (unitized) together with 
ti.j time Interval In which they were taken, are presented below In 
platform coordinates (before realignment). 



■v* 




Sun X. 1 1 

(134:45 - 134:48) +.80/89816 -.00424996 +.58930660 

Moon 

(134:57 - 135:02) -.16146112 -.02178924 -.98663840 


The angle subtended by these two vectors is 135.38 degrees. 


In order to determine the actual subtended angle at the time of 
mark, two sources of data were used. The first was the NASA Apollo 
Trajectory (NAT), from which a spacecraft state vector was obtained 
for an epoch immediately prior to the sightings. This vector and time 
were : 

Time: 17 April 1970, 09:50:00.00 GMT (134:37:00.00 GET) 

X: -.28228044E9 ft 

Y: +.18416553E9 

Z: +.83610975E8 
X: +.533054E4 ft/sec 

Y: - . 581 689E4 

Z: - . 294957E4 

This vector and the TRW Houston Operations Predictor Estimator program 
were used to obtain spacecraft-centered sun and moon vectors from the 
most recent JPL ephemeris tape (JPL 0E69D) at the start and stop times of 
each measurement interval. These vectors were then unitized and averaged. 
They are presented in Basic Reference Coordinates. 

X Y Z 

+.89145748 +.41570609 +.18025559 

-.95321902 +.27595070 +.12338832 

The angle subtended by these vectors is 135.46 degrees. Therefore, the 
actual difference in the angles is seen to be 135.46 - 135.38 « .08 degree. 
Therefore, an LGC calculated angular discrepancy of -1.12 - (.08) * -1.20 
is seen to exist and is attributed to the LGC ephemeris data and software. 


4-20 



Sun 

(134:46.5) 

Moon 

(134:59.0) 







'■ ~~7 

•|V - 


Table 4.2 LM IMU PREDICTED MISALIGNMENTS (DEGREES) 


X 

Y 

Z 

Axis 

Axis 

Axis 


CM IMU Drift 
(At =10 hours) 

-0.200 

-0.040 

-0.020 

LM IMU Drift at Sun Check 
(At = 15 hours) 

-0.045 

-0.022 

+0.382 

Total Drift at Sun Check 
(At « 25 hours; t - 73:47 GET)) 

-0.245 

-0.062 

+0.402 

IMU Drift-Sun Check to TEI 
(At =5.5 hours) 

-0.016 

-0.008 

+0.138 

Total Drift at TEI 

(At r 30.5 hours; t * 79:27) 

-0.261 

-0.070 

+0.540 

Uncertainty Due to CM IMU 
Initial Alignment (la) 

+0.011 

+0.011 

♦0.011 

Uncertainty Due to LM IMU 
Initial Docked Alignment (la) 

♦0.237 

♦0.144 

+0.038 

Total Initialization 
Uncertainty (In) 

+0.237 

♦0.144 

+0.040 

























Table 4.7 DAP CONTROL AXES ATTITUDE ERRORS 






Table 4.B BODY TORQUES CREATED BY TRANSLATIONAL 
POLICIES DEFINED IN GSOP (Reference -1) 


Axis 

Jets 

Body Axis Torque 

T x 

(In-lb) 

T y 

(In-lb) 

T 2 

(In-lb) 

+X (+P) 





4 jet 

2,6,10,14 

0 

-128.2 

+ 1559.2 

2 jet,sys A . 

2,10 

0 

- 64.1 

+ 779.6 

2 jet.sys B 

6,14 

0 

- 64.1 

+ 779.6 

-X (-P) 





4 jet 

1,5,9,13 

0 

-144 

- 1740 

2 jet.sys A 

5,13 

0 

+ 72 

- 870 

2 jet.sys B 

1,9 

0 

+ 72 

- 870 

+Y (+Q) 





2 jet, prl * 

12,16 

+ 72 

0 

-23,874 

1 jet, tack** 

16,15/16,7 

+ 36 ! 

0 

-11,937 


12,3/12,11 

+ 36 • 

0 

-11,937 

! .. -• . 

-Y (-Q) 


j 



2 jet, prl 

4,8 

-72 j 

0 

+23,874 

1 jet, tack 

8,7/8,15 

- 36 

0 

+11,937 


4,3/4,11 

- 36 

0 

+11,937 

+Z (+R) 





2 jet, pri 

7,11 

-870 

+23,874 

0 

1 jet, tack 

7,8/7,16 

-435 

+11,937 

0 


11,4/11,12 

-435 

+11,937 

0 


*Pr1 ■ Primary Mode 
**Tack ■ Tacking Mode 
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Table 4.8 (Continued) 




Body Axis Torque 

Axis 

Jets 

T , 

(in-lb) 

T 

y 

(in-lb) 

T z 

(In-lb) 

-Z (-R) 





2 jet, pri ★ 

3,15 

+870 

-23,874 

0 

1 jet, tack 

15,8/15,16 

+435 

-11,937 

0 


3,4/3,12 

+435 

-11,937 

0 

+U (=+Y , +Z) 





4 jet, pri 

7,11,12,16 

- 798 

+23 ,874 

-23,874 

2 jet, alf* ** 

7,16 

-399 

+11,937 

-11,937 


11,12 

-399 

+11,937 



-11,93/ ! 

-U (=-Y , -Z) 
4 jet, pri 

3,4,8,15 

+798 

-23,874 

1 

+23,874 

2 jet, alt 

8,15 

+399 

-11 ,937 

+11,937 | 


| 3,4 

+399 | 

) 

-11,937 

+11,937 

i 

+V (=-Y, +Z) 
i 4 jet, pri 

i 

i 

4,7,8,11 

-942 

+23,874 

1 

+23,874 

2 jet, alt 

7,8 

-471 

+11,937 

+11,937 


4,11 

-471 

+11,937 

+11,937 

-V (*+Y , -Z) 





4 jet, pri 

3,12,15,16 

+942 

-23,874 

-23,874 

2 jet, alt 

15,16 

+471 

-11,937 

-11,937 


3,12 

+471 

-11,937 

-11,937 


*Pri = Primary Mode 

**Alt * Alternate Mode 
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Table 4.9 BODY TORQUES CREATED BY ROTATIONAL 

POLICIES DEFINED IN GSOP (Reference 4) 

Body Axis Torque 

Jets T^ Ty ^2 

(in-lb) (in-lb) (in-lb) 


I 


4.7.12.15 

4.12 

7.15 
4,7 

7.12 

12.15 

4.15 


+26,440 

+13,220 

+13,2* ** 0 

+12,749 

+12,821 

+13,691 

+13,619 


+11,937 

+11,937 

-11,937 

-11,937 


+11,937 

-11,937 

-11,937 

+11,937 


3.8.11.16 

3.11 

8.16 

8.11 
11,16 
3,16 
3,8 


-26,440 

-13,220 

-13,220 

-13,691 

-13,619 

-12,749 

-12,821 


+11,937 

+11,937 

-11,937 

-11,937 


+11,937 

-11,937 

-11,937 

+11,937 


+ 196.4 + 6,236.8 + 6,187.9 

+ 196.4 - 409.2 + 12.9 

75 



96.4 - 6,229.4 - 6,278.3 

196.4 + 344.6 + 766.7 

0 - 6,574 - 7,045 


* p r \ • Primary Mode 

** Alt » Alternate Mode 



Table 4.9 (Continued) 




Body Axis Torque 

Axis 

Jets 

T X 

t y 

T Z 

1 


(in-lb) 

(in-lb) 

(in-lb) 

+V(= -Y, +Z) 





2 jet, pri 

1,10 

- 166.4 

- 6,229.4 

A 6,187.9 

1 jet-alt 

1 

0 

- 6,574 

+ 6,175 


i 10 

- 166.4 

+ 344.6 

+ 12.9 

i 

1 

- 


— * 

1 -V(»+Y,-Z) 

2 jet. pri 

2,9 

+ 166.4 

+ 6,236.8 

- 6,278.3 

1 jet. alt 

9 

0 

+ 6,646 

- 7,045 

. 

2 

+ 166.4 

- 409.2 

+ 766.7 
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Table 4.10 BODY TORQUES CREATED BY POSSIBLE 
PITCH ROLL ROTATIONAL POLICIES 




Bo<fy Axis Torque 

Axis 

Jets 

T X 

(in-lb) 

t y 

(in-lb) 

T Z 

(in-lb) 

+Y (+Q) 





4 jet 

2,5,9,14 

*■362.8 

+12,473.6 

90.4 

2 jet 

2,5 

+166.4 

+ 6,236.8 

+ 6,941.7 


9,14 

+196.4 

+ 6,236.8 

- 7,032.1 

-Y (-Q) 





4 jet 

1,6,10,13 

-362.8 

-12,458.8 

90.4 

2 jet 

1,6 

-196.4 

- 6,229.4 

+ 6,941.7 


10,13 

-166.4 

- 6,229.4 

- 7,032.1 

+Z (+R) 





4 jet 

1,5,10,14 

+ 30.0 

+ 7.4 

+12,375.8 

2 jet 

1.14 

+196.4 

- 6,983.2 

+ 6,187.9 


5,10 

-166.4 

+ 6,990.6 

+ 6,187.9 

-Z (-R) 


■■ 



4 jet 

2,6,9,13 

B9 

+ 7.4 

-12,556.6 

2 jet 

2,13 

mam 

- 6,983.2 

- 6,278.3 


6,9 

■B 

+ 6,990.6 

- 6,278.3 
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TABLE 4.11 

Sensor Static Bias Errors 


PIC (3/17/70) 

Gyro Static X 0.06°/hr 

Bias Errors Y . 0 .30°/hr 

2 -0.47 tf /hr 


MCC 4 

TEI (4/14/70) (4/15/70) 

*(-0.11 to -0.19)°/hr -0. 12°/hr 

-0.34°/hr -0.23°/hr 

-0.52°/hr -0.61 ’/hr 


Accelerometer X 
Static Bias y 
Errors 

Z 


60.0 ug 
-31.0 ug 

47.0 ug 


-5.0 ng 
-59.0 ug 
57.0 ug 


* 

Reasonable range derived from noisy data 


-16.8 ug 
-59.4 ug 
47.0 ug 


TABLE 4.12 

Sensor Static Bias Stability 




Shift from PIC 
to TEI (30 days) 

Expected 

lc 

Shift from TEI 
to MCC 4 (1 day) 

cxcected 
1 - 

Gyros 

X 

(-0.17 to -0.2S)*°/hr 

0.28°/hr 

(-0.01 to 0.07)* s /hr 

O.LZVnr 


Y 

- 0.04°/hr 

0.23°/hr 

0.11 o /hr 

0. 22 c /hr 


Z 

- 0.05 & /hr 

0.28°/hr 

- 0.09°/hr 

0.2?'', 'hr 

Acceler- 

X 

-65.0 ug 

77 u9 

-11.8 u9 

30 „g 

dneter 

Y 

-28.0 (1 g 

77 u9 

- 0.4 ug 

30 ..g 


Z 

10.0 ug 

77 pg 

-10.0 u9 

30 ■ g 


★ 

Range determined from noisy data. 
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TABLE 4.13 

Equivalent Accelerometer Bias Errors (ug) 



ASA 023 

ASA 023 



Inflight 

Preflight 


Estimate 

Estimate 


(TEI) 

Mean 

3o 

Accelerometer Bias 

X -65 (1) 

0 

135 

and Nonlinearity 

Y -28 (1) 

0 

125 


Z 10 (1) 

0 

147 

X-Scale Factor and 
Dynamic Errors 

X - 1 (2) 

6 

22 

Y and Z Dynamic Errors, 

Y - 4 (2) 

- 3.3 

60 

ASA Accelerometer Internal 
Misalignment and ASA to IM'J 
Mounting Points Misalignment 

Z -58 (2) 

-21.5 

59 

Total (yg) 

X -66 

6 

137 


Y -32 

- 3.3 

138 


Z -48 

-21.5 

158 


NOTES: 1. Inflight Estimate: Shift between freef light period and 

last compensation value. 

2. Inflight Estimate: Difference between measured total 

error and measured fixea bias. 



* * 





TABLE 4.14 

Equivalent Gyro Bias Errors (deg/hr) 





ASA 023 


ASA 023 




Inflight 


Preflight 




Estimate 


Estimate 




(TEI) 









Mean 

3o 

Gyro Fixed Drift 

X 

(- 

.17 to -0.25) 

* (1) 

0 

0.59 


Y 


-0.04 

(D 

0 

0.66 


Z 


-0.05 

(1) 

0 

0.67 

X-Syro Dynamic 

X 

1 

4* 


-0.14 

0.35 

Drift 


(0.0 to 0.08) 

(2) j 



X-Gyro Spin Axis 

V 

A 

i 



0 

0.07 

Hass Unbalance 







Y and Z Gyro Dynamic 

Y 


0.04 

(2) 

-0.08 

0.34 

Dri f t 

Z 


-0.15 

(2) 

-0.09 

0.30 

Total (deg/hr) 

X 


-0.17 


-0.14 

0.69 


Y 


0.0 


-0.08 ‘ 

0.73 


Z 


-0.20 


-0.09 

0.73 

NOTES: 1. Inflight Estimate: 

Shift from last freeflight period 




last compensation value. 


2. Infliqht Estimate: 

Difference 

between measured 

total 




error and measured fixed bias. 

Possible range derived from noisy data. 
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V 




Unfos of Shifts in Dynamic and total Errors to Expectod 1 
aiues of these Shifts at time of TEI 



X 

Y 

Z 

Gyro Dynamic Error 

(1.17 to 1.83)** 

1.18 

-0.60 

Total Gyro Powered Flight 
Error 

-.13 

0.33 

-0.45 

Accelerometer Dynamic Error 
(Y&Z misalignment, 

X scale factor) 

-0.95 

-0.04 

-1.86 

Total Accelerometer 
Powered Flight Error 

-1.58 

-0.62 

-0.50 


* 

★ * 


From ASA 023 Preflight Estimate. 

Range Determined from Noisy Data. Includes spin axis 
mass unbalance. 
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Figure 4-8 TE1 Burn from Prior to 
Bum Initiation through 
Throttle-up-V Axis Phase 
Plane Plot 




Figure 4-9 TE1 Burn Pitch GDA 

Position Through Throttle-Up 
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Mgure 4-12 LM/CSM Conf iguration 

Auto Maneuver and Actitucj 
Hold-U Axis Phase Pianr Mot 





Figure 4-13 LM/CSM Configuration Auto 
Maneuver and Attitude Hold 
V Axis Phase Plane Plot 
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Figure 4-14 LM/CSM Configuration Auto 
Maneuvep-P Axis Phase 
Plane Plot 


























Figure *-)7 IH/CSM Configuration 
Attitude Hold-P Axis 
Phase Plane Plot 












figure 4-19 LM/CSM Configurat ion 
Attitude Hold-V Axis 
Phase Plane Plot 
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AGS/PGNCS Sensed Velocity 
Comparison During TEI Burn 
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Figure 4-29 


Acquisition of Yaw 
Attitude for MCC-4 Burn 
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Figure 4-31 Acquisition of Roll 

Attitude for MCC-4 Burn 
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Figure 4-32 Spacecraft Attitude During 

MCC-4 Burn (AGS Euler Angles) 
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FIGURE 4-34 ACCUMULATED Y-AXIS SENSED VELOCITY 
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79:27:15 28:05 28:55 29:45 30:35 31:25 32:15 

TIME GET (HR:MIN : SEC) 

FIGURE 4-35 ACCUMULATED Z-AXIS SENSED VELOCITY 
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i FIGURE 4-37 Y AXIS VELOCITY DIFFERENCES 

f (NO AGS COMPENSATION) 
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X VELOCITY DIFFERENCE (FT/SEC) 




Y VELOCITY DIFFERENCE (FT/SEC) 
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TIME GET (HR:MIN:SEC) 

FIGURE 4-40 Y AXIS VELOCITY DIFFERENCES (COMPENSATED 
FOR AGS ACCELEROMETER ERRORS) 
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Z VELOCITY DIFFERENCE (FT/SEC) 
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FIGURE 4-41 Z AXIS VELOCITY DIFFERENCES (COMPENSATED 
FOR AGS ACCELEROMETER ERRORS) 
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FIGURE 4-42 AGS/PGNCS ANGULAR DRIFT - X BODY 
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FIGURE 4-43 AGS/PGNCS ANGULAR ‘DRIFT - Y BODY 
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5.0 SEPARATION MANEUVERS 


Prior to re-entry there were two separation maneuvers to extricate 
the CM from the LM/CM/SM stack. 

5.1 SM SEPARATION FROM THE LM/CM 

Service module separation was performed by applying a positive X 
body velocity using the LM RCS thrusters, firing the separation latches, 
and then applying a retro aV using the LM RCS thrusters. The pre-planned 
maneuver called for 0.5 fps of positive aV x followed by a negative 0.5 
fps a V x . LM IMU sensed a V at the time of separation indicated approximately 
0.7 fps was applied in the positive direction and a subsequent 1.9 fps was 
applied in the retro direction. Approximately one minute later the retro 
velocity was decreased by 0.3 fps leaving a net inertial velocity of 
approximately .9 fps applied to the LM/CM stack and 1.6 fps of relative 
velocity between the SM and the LM/CM stack. 

5.2 LM SEPARATION FROM THE CM 

LM separation from the CM was performed by leaving the CM/LM tunnel 
partially pressurized, firing the latches, and allowing the trapped gas to 
force the two vehicles apart. A similar type of separation occurred on 
Apollo 10 between CSM and LM at LM jettison. On Apollo 10, however, the 
tunnel pressure was approximately 80* higher than on Apollo 13. From the 
Apollo 10 data, LM separation velocity was determined and then considering 
the lower tunnel pressure, an impulse was calculated for Apollo 13. 

Applying the analytical impulse to CM and LM, an estimated \V of 2.15 
fps and 0.86 fps was calculated for CM and LM, respectively. Actual experi- 
enced AV at separation was 1.88 fps and 0.65 fps for CM and LM, respectively, 
indicating good agreement with the estimated impulse. The sensed vectors 
at separation are derived below. 



5.2.1 LM AV 


After bias correction, PI PA counts telemetered from the L6C indicate 
the fol 1 owing inertial thrust velocity changes during the 2 second interval 
141:30:00.16 to 141:30:02.16 GET bracketing the time of LM jettison. 

AV X * -.294 ft/sec 

AV y = .555 ft/stc 

AV Z = -.169 ft/sec 

TTie gimbal angles at the time of separation were: CDUX = -156.1707 deg; 
CDUY = -31.4429 deg; CDUZ = 57.4145 deg yielding a direction cosine matrix 
which transforms from platform to body. 



" -0.65 ' 
-0.017 
. -0.003 . 


ft/sec 


AGS telemetry data shows V gx * -.625 during the 1 second interval 
141:29:59.85 GET confirming the LM IMU data. 

5.2.2 CM AV 

Since the CMC telemetry was In low-bit-rate (data available at 1/5 
the normal sample rate) during separation. It was necessary to compute 10 
second DELV's from the available state vector data. Gimbal angles near the 
center of the 10 second computations cycle were used to transform the Iner- 
tial accelerations Into CM Body Coordinates. 

The CM showed evidence of some venting for about 40 seconds after the 
separation event. 
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From the reconstructed thrust accelerations it appears that the CM 
experienced the following velocity change across the separation period: 


a V x (body) * -1.54 ft/sec 
AV y (body) = 0.42 ft/sec 
aV z (body) = 0.99 ft/sec 


The CM DSKY display of thrust velocity in body coordinates (V16 N83) 
during the LM separation maneuver differed significantly from the derived 
values above due to the large Z PIPA bias error that existed. Application 
of this Z-PIPA bias error greatly reduces the discrepancy between the magni- 
tudes of the reconstructed and of the DSKY displayed aV (from 4.19 ft/sec 
to 1.98 ft/sec). The remaining 1.08 ft/sec discrepancy is unexplained. 
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